Are there waves in elastic wave turbulence ? 
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An thin elastic steel plate is excited with a vibrator and its local velocity displays a turbulent-like 
Fourier spectrum. This system is believed to develop elastic wave turbulence. We analyze here 
the motion of the plate with a two-point measurement in order to check, in our real system, a few 
hypotheses required for the Zakharov theory of weak turbulence to apply. We show that the motion 
of the plate is indeed a superposition of bending waves following the theoretical dispersion relation 
of the linear wave equation. The nonlinearities seem to efficiently break the coherence of the waves 
so that no modal structure is observed. Several hypotheses of the weak turbulence theory seem 
to be verified, but nevertheless the theoretical predictions for the wave spectrum are not verified 
experimentally. 

PACS numbers: 46.40.-f,62.30.-|-d,05.45.-a 
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Wave turbulence is a non equilibrium state in which 
dispersive wavetrains interact due to weak nonlinearities. 
The canonical example is the case of gravity waves at the 
surface of the ocean Many other systems belong to 
this class such as capillary waves, nonlinear optical waves, 
superfluids, Alfven or atmospheric waves etc (see for ex- 
ample 0] for a short review). Contrary to the strongly 
nonlinear case of 3D hydrodynamic turbulence, a pertur- 
bative approach to weak turbulence was developed in the 
60-70's. The idea of the theoretical development is, start- 
ing from the nonlinear wave equation, to develop a hierar- 
chy of linear equations for the successive cumulants 
Using the hypothesis of weak nonlinearity, a temporal 
multiscale analysis can be developed which leads to a 
natural closure of this hierarchy. This is the main differ- 
ence from hydrodynamics turbulence for which no such 
rigorous closure exists. From this closure, a kinetic equa- 
tion can be derived describing the slow evolution of the 
Fourier components of the energy spectrum. The energy 
is slowly (compared to the natural frequency of the wave) 
exchanged among modes through the occurrence of res- 
onances between N waves. In the case of gravity waves, 
the resonance involves four waves, and three waves are 
resonant for capillary waves. 

Stationary solutions of these kinetic equations can 
be exhibited Those can be either equilibrium or 

out-of-equilibrium solutions. The former correspond 
to the Rayleigh-Jeans spectrum of equipartition of en- 
ergy among the modes. The typical case where out- 
of-equilibrium solutions are expected is the following: 
the forcing of the system is applied at large scales and 
some additional dissipation mechanism such as diffusiv- 
ity is operating at small scale to dissipate the energy. 
This is also a canonical situation of 3D hydrodynam- 
ics. The kinetic equation conserves energy (and some 
other invariant depending on the case under considera- 
tion). This leads to a direct energy cascade and a self- 
similar Kolmogorov-Zakharov spectrum with a constant 
energy flux across the scales. 

The weak turbulence theory assumes a certain num- 



ber of hypotheses. Among these is the weak nonlinear- 
ity. Another is that the system is supposed to be large 
enough so that no effect of boundaries can be expected 
through the occurrence of cavity modes (the size of the 
system goes to infinity in the derivation) . The goal of this 
article is to check some of these hypotheses in a real sys- 
tem: a thin elastic plate where bending waves are forced. 
Following recent work, this system is believed to be a 
good candidate for wave turbulence We first de- 

scribe our experimental setup and show estimations of 
the spectrum of the normal velocity of the plate. Fi- 
nally, we investigate whether wave structure is indeed 
present in our system, estimate the dispersion relation of 
the observed waves, and relate our observations to some 
hypotheses of the weak turbulence theory. 
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FIG. 1: Experimental setup: a 0.4 mm thick, 1x2 m stainless 
steel plate is hanging under its own weight. It is bolted along 
the full length of one of its short sides on a beamer. Two soft 
springs hold the bottom corners loosely to prevent too large 
excursions of the plate. A vibrator is anchored on the plate 
at about 40 cm from the bottom and excites bending waves 
by moving normally to the plate. The normal velocity of the 
plate is recorded by a laser vibrometer in a region slightly 
above the middle of the plate. 

We report here a study of bending waves of a stain- 
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less steel plate, of size 2x1 and h = 0.4 mm thick, 
held vertically on a short side (Fig. [1]). A electromag- 
netic vibrator type V406 from LDS is fixed at a point 
located 40 cm from the plate bottom. It can move nor- 
mally to the plate to excite bending waves. The typical 
forcing excitation is either sinusoidal (typically at 30 Hz) 
or a low frequency noise (of restricted bandwidth up to 
15 Hz, for example). The forcing is recorded by a force 
probe of type NTC from FGP Sensors and an accelerom- 
eter 4393V from Briiel & Kjaer. The motion of the plate 
is recorded by using a vibrometer type OFV-552 from 
Polytec. The vibrometer is equipped with a two-headed 
sensor that enables us to record the normal velocity of 
the plate at one point or normal velocity differences be- 
tween two points. The measurement points are close to 
the center of the plate or in the upper half of the plate 
so that the motion of the plate is not directly affected by 
the vibrator. 

This kind of thin elastic plate is known to develop 
bending waves following the wave equation [3, 0] 
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where ( is the normal displacement and i? ~ 2 x 10^^ Pa 
is the Young modulus, a ~ 0.3 is the Poisson ratio, 
and p ~ 8000 kg.m^'^ is the specific mass (values are 
typical for stainless steel). A is the Laplacian, and 
{a,b} = dxxa dyyb + dyyU dxxb - 2dxya d^yb. x is the 
stress function and is related to the longitudinal defor- 
mation of the plate Q. The linear dispersion relation 
is 
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For large amplitudes, the wave equation exhibits cubic 
nonlinearities. This nonlinear equation is believed to lead 
to wave turbulence ^ with the following predicted wide- 
band spectrum for the displacement C: 



Ei;{k)=C 
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where P is the energy flux per unit mass (here we use the 
definition of Connaughton et al. 6]), C a dimensionless 
number, and fc* is an ad hoc cutoff wave number. Us- 
ing ^ to change variables from fc to w, the theoretical 
prediction for the velocity power spectrum is then 



E^{uj) = C 
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where uj* — uj[k*) and C is a number. 

Following Connaughton et al. Q, another prediction, 
based on dimensional analysis, can be derived: Starting 



from the linear dispersion relation of the form lo = Afc", 
using the order of the nonlinearity to get the scaling in P 
and assuming power laws, one predicts a velocity power 
spectrum 
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for four- wave interactions (here A = y i2{i^a'^)p ) 

In order to get a more general prediction based on di- 
mensional analysis, one should note that the Young mod- 
ulus appears in the linear part of the wave equation ^ 
combined with h and p but also alone in the nonlinear 
part ([2]). This is different from the case of surface gravity 
waves for example, for which no additional dimensional 
parameters appear in the nonlinear part, due to a differ- 
ent physical nature of the nonlinearity. Here, dimensional 
analysis should use E, alone and not combined in A with 
h and p. In this way, with the parameters P, E, p and 
h (discarding a as without dimension), one predicts that 
the velocity power spectrum should behave like: 
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where g is an unknown function. The theoretical pre- 
diction ([5]) follows the dimensional prediction ^ with 
g{a,b) oc (folna)^/-^. One would then expect the cutoff 
frequency to follow 
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In Q , this high frequency cutoff ui* had to be introduced 
in order to ensure that a non-zero energy flux exists in 
the out-of-equilibrium case. In the dimensional analysis, 
it appears naturally from the nonlinear conservative wave 
equation. Note that in the case of surface gravity waves, 
one has to invoke dissipation (or the crossover to capillary 
waves) to introduce a high-frequency cutoff in the wave 
spectrum. 

In our case, P can be estimated experimentally. The 
mechanical power V transferred to the plate by the vibra- 
tor can be estimated hy V ~ (J'U) , where is the force 
applied to the plate and U the velocity of the plate at the 
vibrator position. For a given plate, P is proportional to 

r. 

Figure [2] displays the velocity power spectrum mea- 
sured at a point near the center of the plate. At the 
highest frequencies, the spectrum decays exponentially 
(thick dashed line). A cutoff frequency can then be es- 
timated by an exponential fit of the spectrum in this 
high-frequency region. As seen in the inset, /* varies as 
P^/'^ in agreement with ([8]). This suggests that the ob- 
served cutoff may not be due to any dissipative process 
but rather to the nonlinear dynamics of the waves (al- 
though some dissipative process eventually absorbs the 
energy flux). 
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FIG. 2: Power spectrum density of the plate normal veloc- 
ity. The frequency is normalized by a cutoff frequency /* 
estimated from the high frequency exponential decay region. 
The spectrum is normalized by V^^'^ . The forcing is sinusoidal 
at 30 Hz and 5 curves are displayed for V — 3.94, 1.94, 1.16, 
0.30 and 0.117 W (from left to right). The dashed line is a 
F~°'^ decay as an eye guide. The thick dashed line is an pure 
exponential decay. Insert: Evolution of the cutoff frequency 
/* with P. The dashed line is a V^^^ fit. 



In the main figure, the frequency has been rescaled 
by /* (i.e., P^^^) and the amplitude by P^/^ showing 
a good collapse as suggested previously by Boudaoud et 
al. The full spectrum, including the high frequency 
cutoff region, seems to follow a master curve. An inertial 
region could be extracted at frequencies just above the 
excitation frequency for which a power law decay can be 
fitted. The best fit exponent is close to -0.6. Altogether 
this leads to a scaling E^iuj) oc po -7^-0.6 .^^jji^jj jg quite 
far from the theoretical prediction Eyioj) oc P^/'^w^ of 
Diiring et al. Q . As discussed by [3| , the scaling exponent 
in P is closer to that of a three-wave interaction. Note 
that the w""'^ scaling is not far from the dimensional 
prediction by Connaughton et al. for a four-wave in- 
teraction. For three waves, their argument would lead to 
w"^/*, so the actual observed exponent lies in between 
the two predictions. Boudaoud et al. suggest that this 
discrepancy of the scaling in P could come from the fact 
that the plate is not ideal, which is certainly true here: 
in particular it is not perfectly flat at rest. This can 
induce quadratic nonlinearities that would cause three- 
wave interactions They also suggest that some sin- 
gular structures like d-cones or ridges may be dominant 
in the dynamics of the plate. 

To tackle the question of the persistence of waves in 
this nonlinear system, we used measurements of the ve- 
locity differences between two points separated spatially 
in the plate, namely 5v(x) = v(x,t) — v{0,t), where x is 
the coordinate along the axis of the rectilinear displace- 
ment of one arm of the two-headed sensor of the vibrom- 
eter (the other arm being at position 0) . Among the hy- 



potheses required for the derivation of the Kolmogorov- 
Zakharov spectrum of weak turbulence, one assumes that 
the wave structure is conserved, that the linear disper- 
sion relation remains valid, and that the size of the sys- 
tem is large enough so that no modal structure due to 
the boundaries is retained. In this spirit, one can assume 
that the motion of the plate is governed by a superpo- 
sition of pure sine waves that remain coherent over the 
length of the measurement so that one can write for the 
Fourier component of Sv at frequency lo: 
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where A[6, to) is the amplitude of the wave traveling with 
an angle 6 to the axis of the measurement. The frequency 
spectrum is then 
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If one assumes that correlations among waves require 
an interaction of more than just two waves (as would be 
expected for the four-wave interaction suggested in this 
case) then one can write, with the additional assumption 
of isotropy: 

{A{e,u)A*{e\oj)) = {\A{oj)f)S{e - 0') . (11) 

Using (fTTj) in (fTO|) . one obtains: 

{\Sv{x,u;)\^) ^ Ank^{\A{uj)\^)il - Mkx)) . (12) 

The prefactor of the Jq Bessel function is simply twice 
the spectrum of the single point velocity (|i;(w)p) so that 



Esyix,u;) = 2E^{lu){1 ~ Mkx)) 
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This result is well known in the study of diffuse fields in 
acoustics and the Jo function that arises here is related 
to the Green function of the wave in free space 0] ■ 

Here, we measure sequentially Sv{x,t) for values of x 
ranging from 1.8 cm to 50 cm. The temporal spectrum of 
each time series is estimated, and all spectra are gathered 
to obtain E^y {uj,x). Examples of these spectra are shown 
in Fig. [3] for various values of the frequency. A Jq Bessel 
function is then fitted to the data. One can see that the 
fit is very good for all frequencies. At high frequencies 
(above 2 kHz typically), the spatial oscillation is not ad- 
equately resolved by our measurement spatial spacing of 
1 cm and sufficient statistical convergence may not be 
attained, as the amplitudes decay quite fast for these fre- 
quencies. 

The spatial scaling factor is shown as a function of the 
frequency and compared to the theoretical wavelength 
in Fig. m The fitted amplitude is also displayed and 
compared to an estimate of the spectrum measured at 
a single point. The estimated dispersion relation is very 
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FIG. 3: Mixed Fourier spectrum Esv{x,lu) of the velocity dif- 
ferences Sv (dots) compared to a fit by a Jo Bessel function 
(continous line) for various values of the frequency. The forc- 
ing is sinusoidal at 30 Hz. 




frequency [Hz] 



FIG. 4: Evolution of the wavelength as a function of the fre- 
quency. Dots: experimental data, continuous line: theoretical 
value from ([3]). Insert: estimation of the spectrum of the nor- 
mal velocity E^iuo) from a direct measurement (average over 
the same positions as the two point measurements) - solid line 
- and extracted from the fit parameter of the mixed spectrum 
E5y(x,uj) - dots. Same dataset as in fig. El 



of the wavelength estimation (not shown). This means 
that no influence of the nonlinearity is observed on this 
quantity, which is one of the requirements for the weak 
turbulence theory to apply. 

The velocity power spectrum estimated from the fit of 
the mixed spectrum Esv{x,Lij) by the Jq Bessel function 
is compared in the insert of Figure d] with a direct 
estimate of Ey{uj). The latter is obtained by averaging 
the spectrum Ey{uj) over the same positions x used to 
estimate Esv{x, uj). One can see that the two estimations 
are coincident at all frequencies. 

From all these observations, one can conclude that the 
motion of the plate is indeed caused by a superposition 
of waves. The wavelength variation as a function of the 
frequency follows quite closely the theoretical dispersion 
relation and no strong nonlinear effect is observed at this 
level. The spectrum of velocity differences follows very 
closely the behavior expected for isotropic waves in an 
infinite medium and no mode structure is evidenced. All 
these facts support the hypotheses required for the the- 
ory of weak turbulence to apply. Nevertheless a mis- 
match is observed between the measured and theoretical 
spectrum. The scaling of the spectrum with the injected 
power is in strong disagreement with the prediction. The 
analytical functional form of the spectrum is also not 
observed. Our two-point measurement support the fact 
that the dynamics of the plate is due to isotropic waves 
and not to more complex dominant structures like ridges 
or d-cones as suggested in 0, HI. The observation of 
waves also supports the usual change of variable from 
wavenumber k to frequency lu through the dispersion re- 
lation. This is the equivalent for wave turbulence of the 
Taylor hypothesis of hydrodynamic turbulence that re- 
lates space scales to time scales. 
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close to the theoretical one (O . A systematic deviation of 
the experimental data is observed at low frequencies. For 
frequencies over 100 Hz, the mismatch is below 10% and 
over 500 Hz it is negligible compared to all uncertainties. 
The mismatch at low frequency can be due to systematic 
errors, as very few oscillations of the mixed spectrum are 
observed, so that the fitting of the Jq function may cause 
systematic errors. At low frequency, the wavelength is 
also comparable with the size of the plate, so that the 
presence of the springs and the fact that the plate is 
held on the vibrator may influence strongly the dynamics 
of the waves. Note that the estimations performed for 
several amplitudes of forcing show no visible variation 
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